Solitons and diffusive modes in the noiseless Burgers equation 

Stability analysis 

Hans C. Fogedby 

* Institute of Physics and Astronomy, University of Aarhus, DK-8000, Aarhus C, Denmark 

and 

NORDITA, Blegdamsvej 17, DK-2100, Copenhagen 0, Denmark 
(February 1, 2008) 

The noiseless Burgers equation in one spatial dimension is analyzed from the point of view of 
a diffusive evolution equation in terms of nonlinear soliton modes and linear diffusive modes. The 
transient evolution of the profile is interpreted as a gas of right hand solitons connected by ramp 
solutions with superposed linear diffusive modes. This picture is supported by a linear stability 
analysis of the soliton mode. The spectrum and phase shift of the diffusive modes are determined. 
In the presence of the soliton the diffusive modes develop a gap in the spectrum and are phase-shifted 
in accordance with Levinson's theorem. The spectrum also exhibits a zero-frequency translation or 
Goldstone mode associated with the broken translational symmetry. 
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I. INTRODUCTION 

This is the first of a series of papers dedicated to an analysis of the Burgers equation in one spatial dimension. To 
set the stage for our main objective, namely the analysis of the noisy Burgers equation for the stochastic growth of an 
interface driven by annealed white noise to be presented in the following paper, we have found it convenient to begin 
with a brief analysis of the deterministic or noiseless Burgers equation. The reason is that many of the features of 
the noiseless case such a nonlinear soliton modes etc. also play a decisive role in the noisy case. The noiseless Burgers 
equation has the form 

— = i^V^u + AmVu . (1.1) 

Here v is a, damping constant controlling the strength of the linear dissipative term. The parameter A characterizes 
the strength of the nonlinear mode coupling term. 

For A = — 1 the equation (1.1) was originally proposed by Burgers |Q| as a model for irrotational or vorticity-free 
hydrodynamics in order to describe one dimensional turbulence . In the present context for general A we consider 
Eq. (1.1) as providing a description of the slope field u = V/i for a growing interface governed by the noiseless or 
deterministic Kardar-Parisi-Zhang equation |^ 

f^=^v'h + \[yhf . (1.2) 

The Burgers equation (1.1) has been the subject of numerous studies as an initial value problem, e.g., ref. 1^, and is 
still of considerable current interest [Q. 

It is an interesting property of the Burgers equation (1.1) that it is exactly soluble in the sense that the nonlinear 
Cole - Hopf transformation 18,0] 



w = exp 
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allows for an exact mapping onto the linear diffusion equation 



'Permanent address 



dw „o 
at 

which can be simply analyzed. For given initial data u{x, t — Q) — uq{x) we thus have 
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where G{x,t) is the Green's function for the diffusion equation (1.5) 



uodx" 



(1.6) 



G(x,t) = [Airvt]-^/^ 
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and u is given by (1.4). 

The Cole - Hopf transformation (1.3-4) thus permits a rather complete analysis of the Burgers equation. The 
relaxational dynamics of the equation is controlled by solitons connected by smooth regions; in the inviscid limit 
V ^ the solitons become shocks connected by ramps |^,^-^. Although the nonlinear character of the equation 
prevents the application of a superposition principle we can still from a qualitative point of view envisage that 
an initial configuration Uf){x) "contains" a certain number of solitons or smoothed shocks. As time progresses the 
configuration passes through a transient regime dominated by a gas of propagating and coalescing solitons with 
superposed linear diffusive modes. At infinite times the configuration eventually decays owing to the inviscid term in 
Eq. (1.1). This qualitative behavior is depicted in Fig. 1. 

It is instructive to compare the nonlinear irreversible and dissipative Burgers equation (1.1) with the nonlinear 
reversible and dispersive evolution equations ||lOt| : 

The equation of motion for the field theory 
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the complex nonlinear Schrodinger equation 



.du 



X\u?u 



and the sine-Gordon equation, 
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= V^u -|- Asinu 



(1.8) 
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(1.10) 



In addition to the linear dispersive modes obtained for A = the above evolution equations all support soliton 
solutions due to the dynamical balance between the linear dispersive term and the nonlinear term controlled by A. 
An initial configuration thus again breaks up into a gas of moving solitons and linear modes. In the case of the 
sine-Gordon and nonlinear Schrodinger equations, the solitons preserve their identity under collisions owing to the 
complete integrability of these systems; this is not the case for the field equation, here the solitons become deformed 
under collisions. 

In the present paper we analyze the Burgers equation (1.1) from the point of view of a soliton-carrying dissipative 
evolution equation. The paper is organized in the following way. In section II wc summarize the general properties 
of the Burgers equation. In section III we discuss in particular the soliton solution and comment on the morphology 
of a growing interface. In section IV we present a linear stability analysis of the Burgers equation and discuss the 
translation mode and the diffusive scattering modes. In section V we summarize our results and present a conclusion. 



II. GENERAL PROPERTIES 



The Burgers equation (1.1) has the form of a nonlinear diffusive evolution equation with a linear diffusive term 
controlled by the damping or viscosity v and a nonlinear mode coupling term characterized by A. In the context of 
fluid motion the nonlinear term gives rise to convection as in the Navier Stokes equation ||l|-^; for an interface the 
term corresponds to a slope dependent growth . 

Under time reversal t —t and the transformation u —> —u the equation is invariant provided v —>■ —v. This 
indicates that the linear diffusive and the nonlinear convective or growth terms play a completely different role. The 
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diffusive term is intrinsically irreversible whereas the growth term, corresponding to a mode coupling, gives rise to a 
cascade in wavenumber space and a genuine transient growth. The transformation t — + — t is absorbed in m — > —u or, 
alternatively, A —A, corresponding to a change of growth direction. 

We also note that the equation is invariant under the parity transformation x — > —x provided u — > —u. This feature 
is related to the presence of a single spatial derivative in the growth term and implies that the equation only supports 
solitons or shocks with one parity, that is right hand solitons. We mention in passing that parity invariance is restored 
in the case of the noisy Burgers equation. This interesting aspect will be considered in the next paper. 

The Burgers equation is also invariant under a more fundamental symmetry, namely the Galilean symmetry group. 
In fact, boosting the equation to a moving frame with velocity Awq, a; — > a; — Aupt it is easily seen by inspection, 
using d/dt —f d/dt + AuqV, that the equation remains invariant provided we shift the amplitude u by uq, i.e., 
u — > u + uq. We note that unlike the and Sine-Gordon equations (1.8) and (1-10) which are invariant under a 
linear Lorentz transformation with no change in the field amplitude and the nonlinear Schrodinger equation (1.9) 
which is invariant under a Galilean transformation accompanied by a space and time dependent phase shift in the 
wave function, the Galilean transformation of the Burgers equation is nonlinear in the sense that the slope field is 
also shifted. Furthermore, the nonlinear coupling strength A enters explicitly in the Galilean symmetry group. 

In the absence of the nonlinear growth term for A = the Burgers equation (1.1) reduces to the linear diffusion 
equation (1.5) supporting linear diffusive modes u oc exp {—icot ± ikx) with quadratic dispersion 

uj = -ivk"^ . (2.1) 

An initial plane wave configuration thus decays with an envelope exp {—vk'^t). More explicitly, defining the Laplace- 
Fourier transform 

uik^uj) — j dxdtexp {iLut — ikx)u{x,t)r]{t) , (2-2) 

where T]{t) is the step function, i.e., r]{t) = 1 for t > 0,ri{t) = for t < 0, and r/(0) = 1/2, and denoting the initial 
slope configuration by uo(fc) = u{k,t = 0), we have the solution 

u{k,uj)^ . 2.3) 
—luj + i/k'' 

displaying a diffusive pole given by Eq. (2.1). For the temporal correlations we obtain in particular 

{u{k, t)u{-k, t')) = {uo{k)uo{-k)) exp [-(t + t')iyk'^] (2.4) 

where (• ■ •) denotes an average over the distribution of initial values. 

On the other hand, in the inviscid limit for vanishing damping ^ — > the Burgers equation (1.1) takes the form 

= XuVu (2.5) 
dt ^ ' 

which has an exact solution given by the implicit equation u — F{x + Xut), where F is an arbitrary profile. Since the 
propagation velocity Xu thus increases with the amplitude u it follows that an initial configuration uq = F{x) breaks 
and that shock waves are generated. From the form of the exact solutions it also follows that the shocks develop 
with right parity, i.e., a positive discontinuity in u. As mentioned above this is consistent with the parity breaking 
properties of Eqs. (1.1) and (2.5). Searching for a static solution of the form u — A + Brj{x — xq) we find by insertion 

u{x) ^ \u+\f]{x - xo) (2.6) 

with arbitrary amplitude a moving shock is then obtained by applying a Galilean boost, i.e., x ^ x — Xuot, 

u ^ u + Uq, yielding shock solutions of Eq. (2.5). It is also easily seen that Eq. (2.5) supports ramp solutions of the 
form 

u{x) = const — — . (2-7) 
At 

The general picture that emerges in the inviscid limit — s- is thus that an initial configuration uq ~ F{x) breaks 
up into a series of right hand shocks connected by ramps; this picture is in fact substantiated by a steepest descent 
analysis of the Cole-Hopf transformation (1.3-4) in the inviscid limit Q. The time evolution is similar to the one 
depicted in Fig. 1 except that the solitons are sharp. Furthermore, in the absence of dissipation the transient regime 
extends to infinite times. In terms of the height field h — J udx, the morphology consists of cusps connected by 
convex parabolic segments |5|JlT[|. 

In the presence of damping or viscosity this picture is not radically changed. The damping leads to an overall 
relaxation of the initial configuration where the energy, based on the hydrodynamical definition J u^dx (kinetic 
energy) is mainly dissipated in the shocks; in other words, the nonlinear mode coupling term gives rise to spatially 
confined hot zones for energy dissipation associated with the solitons. 
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III. THE SOLITON SOLUTION 



Although the Burgers equation (1.1) admits an exact solution by means of the Cole-Hopf transformation (1.3-4), 
we find it useful for our purposes to approach Eq. (1.1) as a diffusive nonlinear evolution equation and following the 
corresponding analysis of Eqs. (1.8-10) search for permanent profile soliton solutions ||l^. Setting u{x, t) = u{x — vt) 
where v is the propagation velocity and using d/dt = — wV Eq. (1.1) can be integrated once 



const 



(3.1) 



Furthermore, imposing the boundary conditions u = it±,VM — * for x — > ±oo, appropriate for a single soliton 
solution, and subtracting Eq. (3.1) for x — > ±oo, we obtain the soliton condition 



2v 

+ u_ = - — 
A 

relating the propagation velocity v of the soliton to the boundary values u± . We note that Eq. 
the fundamental nonlinear Galilean symmetry, being invariant under the transformation v — > 
In terms of the boundary values u± we can also express Eq. (3.1) in the form 

Vm — —{u-^. — u){u — U-) 



(3.2) 

(3.2) is consistent with 
V + Alto, u± — > u± — uq. 



(3.3) 



which implies a positive slope of u between the boundary values u± and > w_, corresponding to a right hand or 
positive parity soliton in accordance with the symmetry property discussed in section II. In the static limit v = Eq. 
(3.2) implies u+ = — m_, that is a symmetric soliton and by quadrature Eq. (3.3) yields the static soliton solution 



uq{x) = tanh [ks{x 



(3.4) 
(3.5) 



We have introduced the characteristic wavenumber kg setting the inverse length scale of the static soliton; xq denotes 
the center of mass position. The width of the soliton is of order l/kg and, unlike the or sine-Gordon soliton related 
to the amplitude u-^-. In the inviscid limit 0, the wave number fc^ — > oo and the soliton approaches the shock 
discontinuity given by Eq. (2.6). 

A moving soliton is obtained by applying a Galilean boost x x — Xugt and shifting the profile by uq. In terms of 
the boundary values u± we thus obtain the propagating soliton solution 



uo{x,t) 



tanh 



^(u+ -u_)(x 



vt — Xq) 



(3.6) 



with velocity v given by Eq. (3.2). In Fig. 2 we have depicted a single soliton solution of the Burgers equation and 
the associated smoothed cusp profile for the height field of a growing interface. 



IV. LINEAR STABILITY ANALYSIS 



In order to investigate the properties of the linear diffusive modes in the presence of the nonlinear soliton mode 
we here perform a linear stability analysis. Since a Galilean transformation allows for a boost to a finite propagation 
velocity it is sufficient to consider the stability of the static soliton uq in Eq. (3.4). Setting u = uq + du we obtain 
by inserting in the Burgers equation (1.1) to linear order in Su the equation of motion for the fluctuations about the 
soliton profile 

dSu 

— ^ iyV'^Su + XuoV5u + X{Vuo)5u . (4.1) 
Absorbing the first order derivative term by means of the transformation []T2| 

and searching for solutions with time dependence Su oc exp (—iujt) we arrive at the linear eigenvalue problem 
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iujSu = j/V Sit — vk^ 



cosh [ks{x — xq)] 



Su . 



(4.3) 



This equation has the same form as the one encountered in the hnear stabihty analysis of the sine-Gordon sohton 
Interpreted as a stationary Schrodinger equation Eq. (4.3) describes a particle with energy iuj and mass 1/2;/ 
in the exactly soluble Bargman potential —2/ cosh^ a:. The spectru m is w ell-known and consists of a single bound 
state for w = and a band of scattering states for uj = —iv{k^ 
and the band of scattering states. 



fc^) Jl3| , p^ . In Fig. 3 we have sketched the potential 



A. The translation mode 

The bound state solution for w = has the form Subs 1/ cosh[ks{x — xq)] and using Eq. (4.2) 

1 



Subs oc 



COSh^ [ks{x — Xq)] ' 



UJ = . 



(4.4) 



This zero frequency mode has a particular significance for soliton-carrying systems. Since Subs oc {du^/dx)Sx it is 
seen that the mode actually corresponds to an infinitesimal translation Sx of the soliton without changing its shape. 
This mode thus restores the broken translational invariance resulting from the choice of a particular center of mass 
position Xq for the soliton - a so-called translation mode. Such Goldstone modes are quite generally associated with 
broken symmetries [p^ , p^p^ . 



B. The diffusive scattering modes 



In a similar way the band of scattering states has the explicit form 

exp (ifca:) k -\- ikgia,'nh[kg(x — Xq)] 



Suss oc 



with dispersion 



cosh [ks{x — Xq)] k — iks 

OJ = -iv{k'^ + kl) . 



(4.5) 



(4.6) 



The modes form a continuum of spatially decaying diffusive states scattering off the soliton. We note that the 
dispersion law (4.6) for the diffusive spectrum compared with Eq. (2.1) for the linear case has developed a gap 
vk^g — X^u^/Av depending on the soliton amplitude. For x ±oo we have Suss oc exp {ikx -\- iS{k)) exp (— fcs|a;|) and 
Suss oc exp {ikx) exp (— fcs|a:;|), respectively, where the plane wave part is phase shifted by the amount 



5{k) = 2tan-i (^) 
k 



(4.7) 



The transmission coefficient is thus given by t(k) = exp(i(5(fc), i.e., |i(fc)|^ = 1, yielding the reflection coefficient 
r{k) = since |?'(fc)p = 1 — |i(A:)l^ = 0. Consequently, the sohton acts as a reflectionless transparent potential on 
the diffusive modes which only experience a phase shift. We also note that the bound state energy is given by the 
pole k = iks in the S-matrix (fc -I- iks)/{k — iks); the residue in Eq. (4.5) yielding the bound state in Eq. (4.4). 
Furthermore, imposing periodic boundary conditions in a system of size N we obtain the density of diffusive modes 



N 1 dSik) 



(4.8) 



2n dk 

For the change of density of states owing to the presence of the soliton , Ap ~ p — N/2Tr, we have, inserting Eq. (4.7), 



Ap 



1 



(4.9) 



We note that J pdk — N ~ 1 in accordance with Levinson's theorem, i.e., the band of diffusive modes is depleted by 
one mode corresponding to the zero frequency translation mode. In Fig. 4 we have depicted the phase shifted-diffusive 
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mode and the translation mode. In Fig. 5 we have shown the diffusive spectrum (4.6) and the pole structure in the 
complex wavenumber plane. 

Finally, regarding the stability of the soliton we note that the linear diffusive mode is damped according to Suss oc 
exp (— r(fc)t) with a damping constant r(fc) — v{k^ + fc^). In the long wavelength limit fc — > 0, r(fc) approaches a 
^ui^/Av, i.e., the gap in the diffusive spectrum. This result is, however, not in conflict with the 



usual argument of a vanishing T{k) for k ^ 0, characteristic of a hydrodynamical diffusive mode |15|. The argument 
generally follows from the conservation of the local slope or velocity (momentum) implied by the structure of the 
Burgers equation (1.1) which can be written as a local conservation law 

S = ("») 

with a current density 

j = -Vu-^u'. (4.11) 

In the long wavelength limit fc — > the conservation law (4.10) usually implies that the damping constant r(fc) — > 0, 
i.e., an infinitely long lived mode. This follows from 

lini MM ^± [ u{x,t)dx = - / Vjdx (4.12) 
fc-^O ot at J J 

which vanishes provided that there are no currents on the boundaries x = zLN/2. However, in the presence of a 
soliton there is an incoming current j{±N/2) = — A/2w^ and the mode decays with a finite life time in the long 

wavelength limit, corresponding to a gap in the spectrum of iuj. 



V. SUMMARY AND CONCLUSION 



In the present paper we have conducted yet a study of the well-documented noiseless or deterministic Burgers 
equation regarding it as a nonlinear diffusive evolution equation. Although the nonlinear Cole-Hopf mapping to a 
linear diffusive equation in principle allows for a rather complete analysis of the equation, we have found it useful 
to emphasize the solitonic aspects, drawing on the parallel with other nonlinear equations such as the cj)'^ and the 
sine-Gordon equations. Note, however, that unlike the 0* or sine-Gordon equations, where the soliton owes its 
stability to a balance between the dispersive effect of the linear term, tending to break up a wave packet construction, 
and the cascade effect in wave number space due to the nonlinear term, stabilizing a particular wave packet form - 
the soliton, the Burgers equation is intrinsically dissipative and an initial configuration will eventually decay due to 
dissipation unless energy is fed into the system. In this regard the Burgers soliton is a dissipative structure in that 
it owes its stability to the energy flux fed by the non-vanishing currents entering at the boundaries. The nonlinear 
term generating an inverse cascade in wave number space thus provides the energy transport to the center of the 
soliton where the energy is dissipated and the soliton owes its stability to the interplay between the linear dissipative 
term and the nonlinear mode coupling term. Nevertheless, it is useful to consider the soliton as the fundamental 
elementary excitation in the Burgers equation determining the nonlinear nonlocal relaxational aspects. As follows 
from a steepest descent analysis in the inviscid limit v ^ 0, an initial configuration evolves into a gas of propagating 
solitons connected by ramps. The present linear stability analysis then shows that the linear modes which for A = 
dominate the relaxational by means of diffusion for X ^ become sub-dominant in the sense that they develop a gap 
in the ioj spectrum. 

In a subsequent paper we consider the Burgers equation driven by spatially uniform stochastic noise rather that 
deterministic currents at the boundaries. We find that the solitonic aspects still determine the physics and that the 
soliton becomes a bona fide elementary excitation in the underlying field theory. 
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FIG. 1. Here we show the transient evolution of the slope field it from an arbitrary initial configuration uq. The transient 
morphology consists of propagating right hand solitons connected by ramps with superposed damped diffusive modes. Both 
the solitons and diffusive modes transport energy which is dissipated predominantly at the soliton positions. At long times the 
profile decays unless it is driven by currents at the boundaries, corresponding to non-vanishing slope. 
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FIG. 2. In a) we show a single moving soliton profile propagating to the left and in b) the corresponding smoothed cusp 
in the growth profile. This configuration is driven by currents at the boundaries, corresponding to non-vanishing u± and is 
persistent in time. 



V, 



FIG. 3. We here sketch the reflectionless Bargman potential fksll — 2/cosh^ [ks{x — xo)]]- We also show the associated 
zero-frequency bound state (BS) or translation mode and the band of scattering states above the gap I'fcf in the spectrum. 
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FIG. 4. In a) we show the phase shifted diffusive mode superposed on the static soliton, indicated by a dashed lino. In b) 
we depict the translation mode giving rise to a solid displacement of the soliton without changing its shape, thus lifting the 
broken translational invariance. 
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K 



FIG. 5. In a) we show the spectrum of diffusive modes in the presence of a soliton. Unlike the quadratic dispersion for the 

linear diffusion equation, the spectrum here exhibits a gap z/fc^ depending on the soliton amplitude. In b) wo show the pole 
structure in the complex k plane. The real ax;is corresponds to the band of scattering states. The pole in the S-matrix on the 
imaginary axis corresponds to the bound state or translation mode. 
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